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We investigate quantum control of an oscillator mode off-resonantly coupled to an ancillary qubit. 

In the strong dispersive regime, we may drive the qubit conditioned on number states of the oscillator, 
which together with displacement operations can achieve universal control of the oscillator. Based 
on our proof of universal control, we provide explicit constructions for arbitrary state preparation 
and arbitrary unitary operation of the oscillator. Moreover, we present an efficient procedure to 
prepare the number state |n) using only O (y/n) operations. We also compare our scheme with known 
quantum control protocols for coupled qubit-oscillator systems. This universal control scheme of 
the oscillator can readily be implemented using superconducting circuits. 
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As an important model for quantum information pro¬ 
cessing, the coupled qubit-oscillator system has been 
actively investigated in various platforms, including 
trapped ions |Tj, nano-photonics [2], cavity QED [3], and 
circuit QED jj. Due to its convenient control, the phys¬ 
ical qubit is usually the primary resource for quantum 
information processing. Meanwhile, the oscillator serves 
as an auxiliary system for quantum state transfer and 
detection [5]. In fact, the oscillator, associated with 
the phononic or photonic mode, may have long coher¬ 
ent times mini and the large Hilbert space associated 
with the oscillator can be used for quantum encoding 
IMS] and autonomous error correction with engineered 
dissipation m- These crucial features call for deeper in¬ 
vestigations into quantum control theory of an oscillator. 

The seminal work by Law and Eberly m has trig¬ 
gered many theoretical and experimental investigations 
to prepare quantum states of the oscillator assisted by 
an ancillary qubit with Jaynes-Cummings (JC) coupling 
P IMIS], while the general problem of implementing 
arbitrary unitary operation remains an outstanding chal¬ 
lenge. Even with recent advances, protocols for universal 
control require either a large number of control opera¬ 
tions m or a more complicated model with an ancil¬ 
lary three-level system El- Meanwhile, development in 
superconducting circuits acting in the strong dispersive 
regime opens new possibilities for universal control of the 
oscillator | T8j . 

In this Letter, we provide schemes for arbitrary state 
preparation and universal control of the oscillator as¬ 
sisted by an ancillary qubit. These schemes utilize the 
dispersive Hamiltonian m along with two types of drives 
associated with the qubit and the oscillator, respectively. 
The key is the capability to drive the qubit im nm [HH2T] 
and impart arbitrary phases conditioned on the number 
state of the oscillator. 

The Hamiltonian of the qubit-oscillator system is 

H = H 0 + H 1 +H 2 , ( 1 ) 
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Figure 1. Energy level diagram of qubit-oscillator system. 
In the rotating frame of the oscillator, the states {|3,u)} n 
have the same energy. After each operation, the population 
(orange circles) remains in the subspace associated with | g). 
The displacement operation (red dashed arrows) couples the 
states | g,n— 1) and \g,n) with strength y/ne for all n. The 
SNAP gate (blue solid arrows) can simultaneously accumu¬ 
late different Berry phases {0 n } to states {\g, n}}. The Berry 
phase 0 n is proportional to the enclosed shaded area in the 
corresponding Bloch sphere, achieved by resonant microwave 
pulses with frequency u> q — n\ (blue oscillatory fields). 

with dispersively coupled qubit and oscillator m 

H 0 =u q \ e)(e | +uj c h - x I e)(e | n, (2) 

time-dependent drive of the oscillator 

Hi = e(t) e lulot a^ + h.c., (3) 

and time-dependent drive of the qubit 

H 2 =n(t)e iu ° t \e){g\+h.c. (4) 

where u> q is the qubit transition frequency between | g) 
and \e),ui c is the oscillator frequency, a^(a) are the raising 
(lowering) operators, h = a^a is the number operator 
of the oscillator, \ is the dispersive coupling, fi(i) and 
e(t) are the time-dependent drives of the qubit and the 
oscillator, respectively. The eigenstates of Hq are \g,n) 
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and |e, n) with oscillator excitation number n = 0,1, • ■ •, 
as illustrated in Fig. [T] 

We consider control schemes with three constraints to 
achieve universal control of the oscillator: 

1. The oscillator and qubit are never driven simulta¬ 
neously (i.e. e(t)Q(t) = 0 for all f); 

2. The qubit is in the ground state |g) whenever the 
oscillator drive is on (i.e. when eft) ^ 0); 

3. The qubit drive is weak compared with the disper¬ 
sive coupling, i.e. |fi(i)| -C 

With the above constraints, we have two types of oper¬ 
ations. The type ® operation (based on Eq. ^ under 
constraint T2) is a displacement operation 

® D[a ) = exp (aa^ — a*a) , (5) 

with a = i f e(t)dt, which can coherently pump or re¬ 
move energy from the oscillator. The type (2) operation 
(based on Eq. Q under constraint ff 3) is a qubit rota¬ 
tion conditional on the number states |18ti22| , which can 
impart number dependent Berry phases. As illustrated 
in Fig. [T] there is a series of transition frequencies of the 
qubit, {uj q — x n } r ,_o! depending on the excitation num¬ 
ber of the oscillator. We can achieve unitary rotations 
between the selected levels {| g,n) , |e,n)} with negligible 
effect to the rest of system, if we drive the qubit with 
fl(t) = n n (t) e~ mxt and |fi n | <C X fT8H2T] . Hence, we 
can impart a Berry phase to a selected number state: 
\g,n) —» e* e " \g,n), with 6 n proportional to the solid an¬ 
gle subtended by the path in the Bloch sphere associ¬ 
ated with {| g,n) , |e, n)} (as illustrated in Fig. [TJ . Since 
the qubit remains in |g) after the operation, we can ef¬ 
fectively obtain a Selective Number-dependent Arbitrary 
Phase (SNAP) operation 

S n (9 n ) = e ie " !"><"!, ( 6 ) 

which imparts phase 9 n to the number state | n). Since 
the excitation number is preserved during the SNAP 
operation (due to constraint ffl), we may drive the 
qubit with multiple frequency components, 12 (t) = 

( t ) e l ^ ulq ~ xn ' t . These will simultaneously accumu¬ 
late different phases 6 n for different number states and 
implement the general SNAP gate 

OO OO 

@ s (e) = n ^ ( 0 «) = E e ’ e " I ™>< n l» ( 7 ) 

n=0 n=0 

where 9 = {0n}J® o is the list of phases. Since 9 n can be 
an arbitrary function of n, the SNAP gate can simulate 
arbitrary non-linear effects that conserve the excitation 
number. For example, if we choose 9 n oc n 2 , the SNAP 
gate effectively induces a Kerr nonlinearity of the oscilla¬ 
tor. With SNAP gates, we just need to consider real dis¬ 
placement, because any complex displacement a = re 11 ^ 
can be decomposed as real displacement and SNAP gates, 
D(a ) = S(9)D(r)S(—9 ), with 9 n = n<j>{ mod r). 


n 

1 

Otn, 2 

&n,3 

1 — F 

0 

- 0.575 

0.682 

- 0.371 

8.3 x 10 -4 

1 

- 0.313 

0.539 

- 0.316 

6.4 x 10“ 4 

2 

- 0.256 

0.441 

- 0.258 

5.2 x 10“ 4 

3 

- 0.222 

0.382 

- 0.223 

4.9 x 10“ 4 

4 

- 0.198 

0.341 

- 0.200 

4.7 x 10“ 4 

5 

- 0.181 

0.312 

- 0.182 

4.6 x 10“ 4 


(a) 


N 

1 - F 

1 

2.7 x 10“ 5 

2 

1.8 x 10“ 5 

3 

1.1 x 10“ 5 

4 

7.2 x 10“ 6 

5 

5.2 x 10“ 6 

6 

4.0 x 10“ 6 


(b) 


Table I. (a) Fidelity and optimized displacements associated 
with state preparation from |i/^init) = | n) to IV’Enai) = \n + 1). 
(b) Fidelity associated with state preparation from |^i n it) = 

|0>to|Vw>-^rElo l">- 


Proof of universality. To show that the operations 
D(a) and S(9) are sufficient for universal control of the 
oscillator, we first identify p = — i(aJ — a) as a genera¬ 
tor of D (a) for real a, and | Q n = J2n'=o l ri, )( ri, |} as 
generators of S(9). Their commutator is 


Jr j. — 1 



= y/n + 1 (| n)(n + 1 | + | n + 1 ){n |), 


which can selectively couple n and n — 1. This gives the 
group commutator 

D(e)R n (e) D(-e)R n (-e) = exp (i J n e 2 + 0(e 3 )) , (8) 

for small real e and for the SNAP gate 
5({e,...e,0,...}) = R n (e) = 

E”'=o e * £ W){n'\ + E^=„+i W){n'\. For any inte- 
r ~ 'i N ~ 1 f - 'i N ~ 1 

ger N > 0, \J n \ an d \ Q n 1 are sufficient 

t 1 n=0 1 J n=0 

to generate the Lie algebra u(N) over the truncated 
number space spanned by {|n) \n < N }, which implies 
universal control of the oscillator [25H25] . 

Explicit construction of target state. Instead of us¬ 
ing infinitesimal evolutions suggested by the universal¬ 
ity proof, we would like to construct a control sequence 
with finite number of steps to prepare the target state 
\if) = ELo °n \ n ) f rom the oscillator ground state |0). 
First, we observe that the target state can be expressed 
as | ip) = 5|?/5), where S = S ({angle (c n )} n ) is the SNAP 
gate and | ip) = E«=o^ \ n ) wit h c„ = \c n \. Hence, we 
only need to consider the new target state j^) with non¬ 
negative c n in the number basis. We then take the strat¬ 
egy to “unroll” the amplitude via a sequence of interme¬ 
diate states j|i/>„)} , with \Jf n ) = (EEi) 5 "' K)) + 

d n | n) and d n = \JY^n’=n ^n’ > which connects the initial 

state \tfo) = |0) and the target state \i/jn) = |^)- The 
key is to perform a rotation U n £ SO{ 2) that acts non- 
trivially in the subspace spanned by {|n),|n + l)}, so 
that U n d n | n) = c n \n) + d n + i1 n + 1) and consequently 
Unli’n) ~ |^n+l) for 71 = 0, 1, - • • , N. 
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Let us first consider the rotation U n that transfers pop¬ 
ulation from |'i/’init) = |ra) to ['(/’final) = |n + l), with an 
efficient implementation 

U n = D(a { ” ] )R n (tt) D{a { / ] )R n (tt) D(a^), (9) 


where R n (tt) = - £”,=0 \n'){n'\ + £~ = 


n'=n -\-1 


»')( 


n \ is a 


SNAP gate with tt phase shift for number states with no 
more than n excitations. To maximize the state prepara¬ 


tion fidelity, F = 


, we obtain optimized 


Unit) 

displacements ( 01 , 02 , 0 : 3 ) as listed in Table |la| More¬ 
over, we can optimize U n for coherently transfer from 
iV’init) = \n) to IV’finai) = sin (0) \n + 1) + cos (0) | n) for 
0 G [0, |] with fidelity better than 0.999. 

Using these building blocks we can first construct a se¬ 
quence of operators U n £ 50(2) acting on the subspaces 
{17T.) , |n + 1)} so that U n |?i) = cos0„ | n) + sin0„ |n + 1) 
with sin0 n = d '\ + ' . The product C/jv-i ■ ■ • U\Uq pro¬ 
vides an initial guess for state preparation. We then 
combine the displacement D{c// , ' > ) from U n with the dis¬ 
placement D(a { 3 n+1) ) from f7 ra +i to reduce the number 
of parameters. Finally, we optimize over all 2 N + 1 dis¬ 
placement parameters with an initial guess based on the 
result from the previous “local” optimizations [28]. Per¬ 
forming the complete procedure requires N optimizations 
over 3 parameters and one optimization over 2N + 1 pa¬ 
rameters, with optimized fidelity F > 0.999 as listed 


in Table lb for state preparation from IV’init) = |0) to 


final) — 


r; =0 \n). 


/N +1 ^n= 0 I 

Sublinear scheme to prepare number state. The above 
scheme of employing 50(2) rotations to unroll the am¬ 
plitudes works generically for arbitrary target states, 
with the number of gates scaling linearly with the high¬ 
est excitation number. However, certain states with 
a narrow distribution of photon numbers can be pre¬ 
pared more efficiently by taking advantage of the fast 
experimentally-available displacement operations. For 
example, the preparation of the number state | n) re¬ 
quires O (n) sequential SO ( 2) rotations from |0) using 
the generic scheme. In contrast, if we start from the co¬ 
herent state D(a) |0) = |a) with a = yfn, which has pop¬ 
ulation distribution centered around |n) with a spread of 
O (y/n), we need only O [y/n) rounds of 50(2) rotations 
to “fold” the coherent state |ct) to the number state | n). 
Fig- ! compares the number of SNAP gates needed be¬ 
tween the generic linear scheme (with O (n) operations) 
and the specialized sublinear schemes (with O {y/n) op¬ 
erations) designed for preparation from | 0 ) to |ro), with 
various target fidelities. For n > 8 , it becomes advanta¬ 
geous to use the specialized sublinear scheme instead of 
the generic scheme. 

Explicit construction of target unitary. A more gen¬ 
eral version of the problem of arbitrary state prepa¬ 
ration is the construction of an arbitrary unitary op¬ 
eration of the oscillator. The first step is to con¬ 
struct the SO{ 2) transformation in the subspace spanned 



Figure 2. The number of SNAP gates for the preparation 
from |0) to | n) (with fixed lower bound of fidelity) for generic 
linear scheme (black diamond line) and specialized sublinear 
schemes with different target fidelities (colored lines). The 
specialized schemes operate on D{y/n ) |0) by folding all pop¬ 
ulation from the subspace {|n — An) ,..., \n + An)} to the 
number state |n). 


by ||n), | n + 1)}, based on the similar construction as 
Eq. ([9]). Different from state preparation discussed ear¬ 
lier, here U n has to fulfill both conditions: U n \n) = 
cos 0 n \n) + sin 9 n \n+ 1 ) and U n \n+1) = — sin 9 n |n) + 
cos 9 n 1 71 + 1) with 50(2) rotation angle 9 n . In addition, 
it is important to impose the constraint ot\ + oli + 0:3 = 0 
to Eq. <§ to minimize the undesired effects to the sub¬ 
space associated with \n’) 7 ^ |n) or \n+ 1). Numerically, 
we optimize the unitary operation fidelity, -Fumtary = 
where N c is the cutoff dimension 


1 

N c 


Tr 


(ytfi n 

,ideal 


(i.e. the size of the matrices used to represent the oper¬ 
ators) |26j|2Z]. The numerical optimum is attained with 
(ai,a 2 ,a 3 ) = (a,-2a,a). 

To construct an arbitrary unitary Utarget in the 
{ 10 ),..., |tx — 1 )} subspace, we start by taking its inverse: 


U t 


-1 

target 



where W n is the non-trivial block and J n < is the n! x n! 
identity matrix. We first apply a SNAP gate such that 
the last column of the W n block now contains only non¬ 
negative coefficients. We then apply n — 1 consecutive 
50(2) rotations eliminating the off-diagonal elements in 
the last column of W n , such that the column becomes 
(0, • • • ,0,1)^. Since all rows of a unitary matrix are or¬ 
thonormal, the last row of the W n block must be trans¬ 
formed to (0, • • • , 0, l) n . Hence, the result is 


T/( n ) f/( n ) T/( n ) Q(n)fr—T _ 

v n -1 v n-2 ‘ ‘ ‘ v 0 J ^target — 


(Wn-! 

0 

0 ) 

0 

1 

0 

V 0 

0 

IN c — nj 


where 5 ^ is a SNAP gate necessary for any complex 
phases unobtainable with the 50(2) operations. The n — 
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(a) 




(C) 


0 2 4 6 



(d) 



Figure 3. (a) An example target unitary operation, [/target 
(a permutation). (b) The product is close to identity, 
[/construct ^target ~ T with small non-zero off-diagonal ele¬ 
ments. (c,d) Comparison of fidelities after column-wise op¬ 
timization (black) and fidelities with additional optimization 
over all displacement parameters (red) for (c) Fourier (trian¬ 
gles) and permutation (squares) operations and (d) randomly 
generated target unitary operations (hexagons). 


1 SO(2) rotations 


Vk n> = D{ar)Rk M D(-2a); l ’)Rk (tt) !>«') 

can be individually optimized, before being chained to¬ 
gether for a second round of optimization over n — 1 dis- 

f / ii n ~ 1 

placement parameters { ax' > . The cost function to 

l J k=o 

be minimized for the second round of optimization is the 
sum of absolute values of off-diagonal terms (excluding 
W n -1 block). We iterate the procedure until we obtain 


(™)\ 


„( n )\ 


U t 


construct ^target 


( 10 ) 


With ^construct = 1^=1 ‘ ^ ^^ 

illustrated in Fig. I a ,b) for a specific U tal get [28, ■ 

Using the above decomposition, we need n[n — l)/2 
SO( 2) rotations (each containing 3 displacements and 
2 SNAP gates). We can combine consecutive displace¬ 
ments, lowering the number of displacements to 2 per 
SO{ 2) rotation. We also need one SNAP gate at each 
iteration of the W n -A W n -i step, with a total of n addi¬ 
tional SNAP gates. For various [/target, as illustrated in 
Fig.§c ,d), we find that the step-wise optimization pro¬ 
cedure can yield good final fidelity F un i tar y > 0.99, which 
can be further improved to -F U nitar y > 0.999 with a third 
round of simultaneous optimization over all n(n — 1 ) /2 

displacement parameters {a). }k<n'<n- 

Our scheme can be applied to the general Hamiltonian 
with dispersive coupling, — % n |e, n) (e, n |, as long as 
the number-dependent qubit frequency shift can be re¬ 
solved (|Xn ~ Xn'^n\ (t )| , 7 , k, for all relevant n and 


n! , where 7 and k are the qubit and cavity decoherence 
rates, respectively). Furthermore, we can also extend 
the arbitrary unitary control to the subspaces spanned 
by {| g,n!) , |e, n , )}™ /= 1 0 , so that we can control the entire 
2 n-dimensional quantum system. 

Discussions. We now compare our SNAP-gate-based 
quantum control scheme with previous protocols. The 
scheme proposed by Law and Eberly [12] is based on the 
JC model, iTjc oc aa + + h.c., which enables preparation 
of arbitrary superposition of number states. The scheme 
by Mischuck and Molmer m further extended the JC 
model from state preparation to arbitrary unitary opera¬ 
tion, but it is rather complicated because any JC control 
pulse necessarily couples the states \g,n) and \e,n— 1 ) 
(for all n ) simultaneously and with varying strength g^/n. 
In contrast, our scheme is based on the dispersive qubit- 
oscillator coupling, iJdispersive = ~X l e ) ( e l which pre¬ 
serves the oscillator number states, enables the SNAP 
gate to directly access the two selected sublevels \g,n) 
and |e,n) with negligible effects to the rest of the levels, 
and ultimately leads to efficient universal control of the 
oscillator. Similarly, the proposal by Santos m intro¬ 
duces a different model with a three-level A-type ancil¬ 
lary system to achieve universal control, but it is experi¬ 
mentally more challenging than the simple qubit ancilla 
considered in our scheme. 

With dispersive qubit-oscillator coupling, there are 
other control protocols available. For example, in the 
presence of the oscillator drive e, we may “block” the 
processes \n' ± 1 ) — > \n’) by driving the qubit resonantly 
to the transition | g,n') -f-> \e,n'), with ey/W 

|2~9] . Similarly, by resonantly driving transitions | g, n') <A 
| e,n’) for n! < n and n! > n + 1 , we block all number 
changing transitions, except for the transition between 
{|n), |n + 1)} that can be used for SO (2) unitary rota¬ 
tions. This blockade scheme is relatively slow and each 
elementary operation takes time r ~ (ey/n) fl _1 
X -1 due to the blockade requirements, while the SNAP- 
gate-based scheme can be much faster with r ~ x -1 
by relaxing contraint ^3 and numerically optimizing the 
shaped pulses with fl (t) ~ X- 

In conclusion, the SNAP-gate-based scheme provides 
universal control of the oscillator mode with strong dis¬ 
persive coupling to a qubit. Based on the proof of uni¬ 
versal control, we show explicit constructions for arbi¬ 
trary state preparation and arbitrary unitary operation 
of the oscillator. We also present an efficient procedure 
to prepare the number state | n) using only 0{y/n ) op¬ 
erations. We note that deterministic SNAP-gate-based 
preparation of \n = 1 ) photon number state has been 
demonstrated using superconducting circuits [30!. The 
techniques introduced here are not restricted to oscil¬ 
lator modes such as mechanical motions [T] and opti¬ 
cal/ microwave cavities EH5, but can be extended to 
multi-level systems such as Rydberg atoms with large 
angular momentum [29., as long as the dispersive cou- 
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pling between the qubit and oscillator/multi-level system 
is strong. 
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